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Abstract 

A stationary subdivision sciieme generates the full space of polynomi- 
als of degree up to k if and only if its mask satisfies sum rules of order 
k + 1, or its symbol satisfies zero conditions of order fc + 1. This property 
is often called the polynomial reproduction property of the subdivision 
scheme. It is a well-known fact that this property is, in general, only nec- 
essary for the associated refinable function to have approximation order 
A: + 1. 

In this paper we study a different polynomial reproduction property of 
multivariate scalar subdivision scheme with dilation matrix ml, \m\ > 2. 
Namely, we are interested in capability of a subdivision scheme to repro- 
duce in the limit exactly the same polynomials from which the data is 
sampled. The motivation for this paper are the results in [7] that state 
that such a reproduction property of degree k of the subdivision scheme 
is sufficient for having approximation order k + 1. 

Our main result yields simple algebraic conditions on the subdivision sym- 
bol for computing the exact degree of such polynomial reproduction and 
also for determining the associated parametrization. The parametriza- 
tion determines the grid points to which the newly computed values are 
attached at each subdivision iteration to ensure the higher degree of poly- 
nomial reproduction. We illustrate our results with several examples. 

Keyword: Subdivision schemes, polynomial reproduction, subdivision parametriza- 
tion, approximation order 



Introduction 



Interest in subdivision schemes is motivated by their applications in computer 
graphics, computer aided geometric design, animation, wavelet and frame con- 
struction. Important properties of subdivision schemes such as convergence, 
regularity, polynomial generation, approximation order, etc., have been stud- 
ied by several authors, see surveys [3l HI |8l IH [TH [15] and references therein. 
In this paper we would like to distinguish between the concepts of polynomial 
generation and polynomial reproduction of subdivision schemes. The so-called 
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polynomial generation of degree k is the capability of subdivision to generate 
the full space of polynomials of degree up to k. This property is equivalent to 
sum rules of order fc + 1 on the subdivision mask, or, equivalently, to zero con- 
ditions of order 1 on the subdivision symbol, see e.g. [IKHIITI]. Polynomial 
generation of degree k also implies that the associated refinable function has 
accuracy of order k [TS], but is, in general, only necessary for the corresponding 
shift-invariant space to have approximation order k + I. This has been already 
observed for B-splines and box splines [T]. The so-called polynomial reproduc- 
tion is the capability of subdivision schemes to produce in the limit exactly the 
same polynomials from which the data is sampled. The results in [71 Section 
2.4] state that polynomial reproduction of degree k of convergent subdivision 
is sufficient for the associated shift-invariant space to have approximation order 
fc -t- 1. This motivates our interest in polynomial reproduction of subdivision 
schemes. 

Our main goal is to derive simple algebraic conditions on subdivision symbol 
that allow us to determine the degree of its polynomial reproduction. Note that 
the concepts of polynomial reproduction and generation coincide in the case 
fc = 0, i.e. in the case of reproduction of constants. In the L2 setting, poly- 
nomial generation of degree fc is also sufhcient for approximation order fc -I- 1, 
see e.g. [12]. For convergent interpolatory subdivision schemes, the concepts of 
polynomial generation and reproduction are equivalent and, thus, characterize 
the approximation power of the corresponding shift- invariant space, see 
We emphasize that there is a multitude of results on approximation order of a 
refinable function. Those results however are mostly derived from the proper- 
ties of the associated shift-invariant space and not from the properties of the 
coefficients of the refinement equation - the subdivision mask, in the subdivision 
context. There are also data pre-processing techniques for achieving the opti- 
mal approximation order of a shift-invariant space associate with a convergent 
subdivision scheme, [3 |9l [13] . The polynomial reproduction of order fc makes 
preprocessing unnecessary, which is undoubtedly an advantage. 

In the univariate case, polynomial reproduction has been studied in |10] for 
binary primal and dual subdivision schemes and extended in [S] to univariate 
subdivision schemes of any a-rity. In [3] the authors provide unified algebraic 
conditions on the subdivision symbol for polynomial reproduction with no re- 
strictions on the associated parametrization - the grid points to which the newly 
computed values are attached at each subdivision iteration. The results of our 
paper extend [6] to the multivariate setting for scalar subdivision with dila- 
tion matrix ml, \m\ > 2. To the best of our knowledge our results are the 
first ones on multivariate polynomial reproduction of subdivision schemes. Our 
interest in the case of dilation matrix ml is motivated, e.g., by the bivariate 
•\/3— subdivision whose refinable function is also refinable with respect to dila- 
tion —3/ and iterated mask 0(2:122'^, z\z2^) ■ a{zi, Z2). There are several other 
expansive dilation matrices M which satisfy M" = m/, and, thus such that our 
results are applicable. 

The main result of our paper. Theorem 12. 6[ states that for a non-singular 
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subdivision scheme with finitely supported mask a = {ao,, a G Z*} and symbol 
a{z) — aa z" the polynomial reproduction of order k is equivalent to 

(£)Ja)(l) = |mr n n - and = for e e S', |j| < fc , 

(0.1) 

where S' is a finite set of certain multi- indices and t = (ti,--- ,Ts) G K** 
appears in the parametrization associated with the subdivision scheme. The 
importance of condition (jO.ip for fc = 1 is that it allows us to identify the 
correct parametrization that for any non-singular or for even only convergent 
subdivision scheme guarantees at least the reproduction of linear polynomial. 
The parametrization determines the grid points to which the newly computed 
values are attached at each step of subdivision recursion to ensure the higher 
degree of polynomial reproduction of a scheme. 

This paper is organized as follows. The first section sets the notation, pro- 
vides the background on multivariate subdivision schemes stressing the differ- 
ence between polynomial reproduction and polynomial generation. We also 
define sequence of parameter values associated with a subdivision scheme, i.e., 
the subdivision parametrization. In section 2, we first provide algebraic tools for 
determining the correct parametrization needed to ensure reproduction of linear 
polynomials. We also give the necessary and sufficient conditions on the symbol 
of any non-singular subdivision scheme that guarantee polynomial reproduction 
up to a certain degree. Effect of a shift of the mask on the degree of polynomial 
reproduction is investigated in Section 3. There we also show that the concept 
of polynomial reproduction and polynomial generation are equivalent for con- 
vergent interpolatory schemes. Thus, reproducing the results in [11' . In section 
3 we also provide the correct parametrization for box spline subdivision schemes 
together with several examples. The effect of the shifts of the box splines on 
the approximation order of the corresponding shift-invariant spaces has been 
already observed e.g. in [T]. 

1 Background and notation 

1.1 Subdivision scheme 

A scalar s-variate subdivision scheme with a dilation matrix ml, |to| > 2, is 
given by a scalar finitely supported sequence a = {ao, G R, a G Z*^}, the so- 
called mask. The subdivision operator iSa acting on data sequences d = {da, a £ 
Z"} G £(Z^) is defined by 

(5ad)^ = ^ a.a-,np d^ , a G Z'* , (1.1) 
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where is the space of scalar sequences mdexed by Z'*. A subdivision scheme 
is the recursive algorithm given by 

d('-+i) = ^ad^"^) = J2 a«-m/3dj^ , d(") e £{Z'), re No, (1.2) 

where Nq is the set of natural numbers including zero. 

The symbol of a subdivision scheme is given by the Laurent polynomial 

a{z)^ a„^°, z = {zi,...,z,)e{C\{0}y , 

where for cx = (ai , . . . , a^) G Z'* we define z" ~ z"^ Z2 ^ • . . . • zf" . Denoting with 

E = {0,...,|m|-ir, (1.3) 

the set of representatives of 1/ jml/ containing = (0, 0, ... , 0), the \m\'^ sub- 
masks and their symbols a(,{z) are defined by 

{ae+™„, a e Z^} (1.4) 

and 

a,{z) - J2 ' e ^ E , (1.5) 

respectively. Then, we get the following decomposition of the mask symbol 

a(^) = ^ ae(2;). (1.6) 

eeE 

1.2 Sum rules and zero conditions 

The sum rules of order 1 in terms of submasks read as follows 

ae(l) = ^ ae+ma = 1 , 6 G E . (1.7) 

For I = (Ci, . . . , is) e set = e'^''^^ , j = 1, . . . , s, the set in corre- 
sponds to 

S - Se = {£ = ™ : eeE}, (1.8) 
and contains 1 = (1,1,. ..,1). The sum rules of order 1 take an equivalent form 

a{l) = \m\' and a(£) = for e e S' := S \ {1} . (1.9) 

Following the notation in fF, we call S' the zero set, and the conditions in 
(jl.9p the zero condition of order one (Condition Z\). In the literature, both 
the conditions in (|1.7p and their equivalent form in (|1.9|) are called the sum 
rules of order one. We also make use of the higher order sum rules, see [12] and 
references therein: The mask symbol a.{z) is said to satisfy the zero condition 
of order k ( Condition Zk ), if 

a(l) = |m|" and {iya){e) = for £ e S' and |j| < fc . (1.10) 

We denote by the j-th directional directional derivative. 
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1.3 Parametrization 



Since most of the properties of a subdivision scheme, e.g. its convergence, 
smoothness or its support size, do not depend on the choice of the associated 

(r) (r) 

parameter values , a G , to which the data d^ , a G Z'* , generated by the 
r-th step of a subdivision recursion is attached, these are usually set to 

tW:=-^, aeZ^ r>0. (1.11) 

We refer to the choice in (jl.lip as standard parametrization. We show in section 
[5] that the capability of subdivision to reproduce polynomials does depend on 
the choice of the associated parameter values and the standard parametrization 
is not always the optimal one. As in [B^ for the univariate case, the choice 

t(:':=t|,^^ + 4' t(^'=tr''-^- C=0. "er, r>0, (1.12) 

nv 

with a suitable r € R'' turns out to be a better selection. 

We call the sequence {t''"\ k > 0}, with t'*"^ = {ta\ a £ V} the sequence of 
parameter values associated with the subdivision scheme. 

1.4 Convergence and non-singularity of subdivision 

Following jlOj . our definition of convergence depends on the parameter values 
associated with a given subdivision scheme. Since the subdivision process gener- 
ates denser and denser sequences of data d^''^ r > 0, a notion of convergence can 
be established by using a sequence {F'''\r > 0} of continuous functions F'^^'^ 
that interpolate the data d'-''-' at the parameter values {t'^'"', fc > 0} associated 
to the subdivision scheme, namely 

F('-)(t([)) =d([\ aeZ^ r>0. (1.13) 

Definition 1.1. If the sequence of continuous functions {F'^'^\ r > 0} satisfying 
(|1.13p converges, then we denote its limit by 

3d := lim F^'^\ 

r— >oo 

and say that g^ is the limit function of the associated subdivision scheme 

for the initial data d*-*^-* = {da\ a. G Z'*} G The limit function (f) '■— gs 

with the initial data 

^" \^ 0, otherwise ' & i 

is called the basic limit function of this scheme, it is compactly supported and 
satisfies the refinement equation 

(j> = ^ aa(l){m ■ —ot) (1-14) 
with refinement coefficients given by the mask a. 
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Definition 1.2. A subdivision scheme is called non-singular, if it is convergent, 
and gd = if and only if d is the zero sequence, i.e. do, = for all ol . 

Next we show that the notion of non-singular subdivision scheme is equiva- 
lent to the notion of linear independence of the integer shifts of its basic limit 
function. 

Proposition 1.3. A convergent subdivision scheme S'a is non-singular if and 
only if the integer translates of the solution <j) of refinement equation \1.14^ with 
coefficients satisfying ae(l) = 1, e g E, are linearly independent and form a 
partition of unity. 

Proof: 

"=>:" Assume that S'a is convergent, then the basic limit function (j) = S^d 
satisfies the refinement equation (jl.l4l) and its integer shifts form a partition 
of unity. If the convergent subdivision scheme S'a is non-singular, then for any 
starting sequence d € £(Z*) we have 

Srd = J2 d«<^(- - a) = 
if and only if d is the zero sequence. 

"-4=:" Assume that the integer translates of the solution of refinement equa- 
tion (jl.l4p are linearly independent and form a partition of unity. This implies 
that the subdivision scheme S'a associated with the symbol a{z) is convergent, 
see [SI Lemma 2.3]. It is also then non-singular, otherwise one easily gets a 
contradiction to the assumption on linear independence of the translates of (j). 



1.5 Polynomial generation versus polynomial reproduc- 
tion 

We denote by life the space of multivariate polynomial of total degree fc e Nq. 

Definition 1.4 (Polynomial generation). A convergent stationary subdivision 
scheme Sa generates polynomials up to degree da (is TldQ- generating) if for 
any polynomial tt G 11^^ there exists soTfic initial ddta, ' G ^oo 

such 

that S^<^'^^ — IT. Moreover, the initial data q*^"^ is sampled from a polynomial 
TT € Ilrfg with the same leading coefficients as n G 11^,^ . 

Note that the assumptions on the properties of tt in the above Definition 
are justified by [71 Lemma 2.1]. Note also that polynomial generation is also 
studied in [J]. 

We continue with a slightly different notion, the notion of polynomial repro- 
duction which requires a specific choice of starting sequences of a polynomial 
limit. The concepts of polynomial reproduction and generation coincide for 
d-G = dR = 0. 
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Definition 1.5 (Polynomial reproduction). A convergent subdivision scheme 
Sa, with parameter values {t'^^\ r > 0} is reproducing polynomials up to degree 
du (is Hdn-reproducing) if for any polynomial tt G ^cIr and for the initial data 
pC) = {7r(ta-'), a £ Z*^} the limit of the subdivision satisfies S'^p'^"-' = vr. 

Another important property of subdivision is the so-called step-wise poly- 
nomial reproduction, we make use of it in Section 2. 

Definition 1.6 (Step-wise polynomial reproduction). A convergent subdivision 
scheme 5a with parameter values {t^'^\ r > 0} is step- wise polynomial repro- 
ducing up to degree k (is step-wise Ilk -reproducing) 'if for any polynomial tt G 11^ 
and for the data d^'') = {7r(tL''^), a e Z'} 

^(r+i) ^ or, equivalently, 7r(t(;+i)) = ^ a„-m/3 7r(t^''^), a e Z^ 

(1.15) 

The next proposition shows that for a non-singular subdivision scheme the 
concepts of polynomial reproduction and step- wise polynomial reproduction are 
equivalent. 

Proposition 1.7. A non-singular, subdivision scheme is step-wise polyno- 
mial reproducing up to degree k if and only if it is polynomial reproducing up 
to degree k. 

Proof: "=4>:" For any polynomial tt e Uk, let d(") := {7r(tL°^), a e Z"}. 
If the subdivision scheme S'a is step-wise n^-reproducing, then the sequence 
{F''^\ r > 0} of continuous functions F'""' satisfying (|1.13l) with da'' = tt (^^a^^ , 
a S Z", converges uniformly to tt as r — ^ oo, since the distance between the grid 
points m~'^Z^ goes to zero. 

"■4=:" Let us assume next that the subdivision scheme 5a is Uk -reproducing. 
Let r > 0. On the one hand, applying the subdivision scheme to the data 
d('') = {7r(tL''^), a e Z"} we obtain, 

On the other hand, for the sequence p'-'"+^' := {7r(ta^^-'), a e Z^} we also have 
Therefore, by the linearity of the operator Sa it follows 

and, thus, p(''+^) — d''"+^^ — 5ad''"^ due to the assumption of non-singularity. 
Thus, the claim follows. ■ 



7 



2 Algebraic condition for polynomial generation 
and reproduction 

In this section, for a non-singular subdivision scheme, we determine the value of 
T € in (|1.12l) that guarantees the polynomial reproduction of linear polyno- 
mials, see Proposition l2.3l In Theorem l2.6[ we then provide algebraic conditions 
on a{z) for checking the reproduction of polynomials of higher degree. These 
algebraic conditions depend on the previously obtained value of r. In the case 
of only convergent schemes see Corollaries 12 .41 and 12 . 71 We start by defining the 
tensor product polynomial of degree |j|, j € Ng, given by 

s J, -1 

qo{z):^l, qj{zi,...,Zs) ■.= Y[Y[{zi- e,), j = (ji, . . . , Js)- (2.1) 

i=l £i=0 

The following auxiliary proposition states known results on polynomial genera- 
tion that we make use of in this section. 

Proposition 2.1. Let S'a be a convergent subdivision scheme. 

(i) The subdivision scheme S'a reproduces constant sequences or, equivalently, 

its symbol a{z) satisfies \1. 7p , ij and only if a(z) satisfies (jl.9p . 

(ii) Let fc G N. The symbol a{z) satisfies condition if and only if 

{Lya^){l) = \m\-'L^a{l), for e G E, j G N§, |j| < fc. 

{Hi) Let /c G N. The symbol a{z) satisfies Condition if and only if 

^ gj(Q;-m/3)a„-™/3 = |m|""Dja(l), a G Z", j G Ng, |j| < /c. 

(2.2) 

Proof: The proof of (z) for dilation matrix 2/ follows from the definition of the 
operator S'a and [4, Section 6] and in the case of dilation matrix to/, |m| > 2, 
is in [13l Lemma 3.3 ]. The proof of {ii): Let j G Ng, |j| < k. Since af,{z) = 
ae-mct^;*^^™", its j-th derivative satisfies 

(i:^ae)(;z) - ^ gj(e-TOa)ae-,„„^^-™°-J, ;s G (C \ {0})^ (2.3) 
Next, due to a{z) = ^ ae(^), we have 

{D^a){e) - ^(i:^ae)(£) = ^ E "^J^^ " "ia)ae-™«£"-'""-j 
= E e"^^ E gj(e - ■mcx)&e_rna = ^ e"^^ {l^ae){l) 

eSE aSZ" eSE 
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for £ G ^. The rows of the matrix of this hnear systera are given by (e^ ■')eeE- 
This matrix is invertible and Condition Z^, due to 



£ = 

eeE 



|to|", £ = 1, 
0, otherwise, 



is equivalent to 

eeE 

Proof of (iii): For any o; G there exists eeE and /3 G such that 
a = e + TO/3. Thus, 

D^af,(X) = ^ qj(a - m/3)aa_m/3, eeE. 
/3ez= 

The claim follows by (m). ■ 

Remark 2.2. A^ote t/iat i/ie proof of (ii) is also implied, for example, by 
Theorem 3. 7]. 

The next result provides a simple algebraic conditions for determining t g 
R*, which appears in ()1.12p and guarantees the reproduction of linear polyno- 
mials. 

Proposition 2.3. Let be a non-singular subdivision scheme that generates 
linear polynomials, i.e. its symbol satisfies Condition Z\. Then S^, reproduces 
linear polynomials if and only if its parameter values are given by hi. 12]) with 
T= \m\-^ {D^^a{l),...,D^-a{l)). 

Proof: According to Proposition ll . 71 for non-singular subdivision schemes, poly- 
nomial reproduction is equivalent to step-wise polynomial reproduction. More- 
over, any convergent subdivision scheme reproduces the constants, hence it is 
sufficient to prove the claim for polynomials of the form -k^xi, . . . ,Xs) = Xj, 

j = 1, . . . , s. Let r e No and set da"* = 7r(ta'), a G Z-*, with Tr{x) — Xj. Then 
for any a e Z* and e = (ei, . . . , e^) we get 

,(r+l) _ ,(r) _ ,(r) _ f (r) «3 " /^J A 

t^/^Ws /q^zT's r^iz^s \ / 




pel 

pel 



Oe(l) 

(r) maj + ej\ D'^^a(,{l) 



TO 



r+1 / m'"+l 



e,\ D-^ajl) 

I™ I s . ™r4-l 
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where the last equaUty is due to Proposition 12.11 part (m) for j = ej. Thus, 
dma+e IS equal to 

^(fir+i) \ _ .('■+1) , "^"J + _ fir) , WOj - Tj + Cj 

T^yi^rna+e) " + ^r+1 " ^O.j + ' " t ^ , 

if and only if Tj — \m\^'^D'^^ a{l). ■ 

Note that not all convergent subdivision schemes are non-singular, e.g. the 
subdivision scheme based on the four directional box spline symbol is not. For 
such schemes one can still determine a parametrization that ensures its poly- 
nomial generation property. The following result is a direct consequence of 
Proposition [ 



Corollary 2.4. Let be a convergent subdivision scheme that generates linear 
polynomials, i.e. its symbol satisfies Condition Z\. Then reproduces linear 
polynomials if its parameter values are given by with 

T^\mr {D^'a{l),...,D^'^a{l)). 

The following result is crucial for the proof of the main result, Theorem l2.6l 
of this section. It allows us to express the polynomial generation of S'a in terms 
of the properties of its symbol. 

Proposition 2.5. Let k E N, t eW and as in (|2.ip . A subdivision symbol 
a{z) satisfies 

(l?Ja)(l) = |m|^qj(T), {iya){e)=0 e G S', j e Ng, |j| < fc, (2.4) 
if and only if 

5] a„-™/3/3J = l"^^]', aeZ^ jeNg, |j|<fc. (2.5) 



Proof: First, note that due to Proposition 12 . II the conditions in (12. 4p are equiv- 
alent to 

= " TO/3)a„_™^, j e Ng, |j| < fc . (2.6) 

Using this equivalent formulation, we prove the proposition by induction on k. 
For k = the claim is true since for any r S M'' we get 

loi^) = J2 = ( ^^-^ ) = 1- 

Next, we assume that the claim is true for all |j| < fc — 1 and prove it for j € Ng 
with |j| = fc. To this purpose, we write the polynomial in x of (total) degree 
|j| = fc as 

qj{a-mx)= ^ cia,£X^, xeW, Cj,a,j7^0. (2.7) 
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Therefore, using the induction assumption and by (12.61) and (|2.7p we have 



^eNg, |£|=fc /JeZ" ^GNg, |£|<fc-l 



£eNg, |£|=fc \/3ez= ^ ^ 

The last equaUty is due to the fact that 

'7j(-r) =gj fa-™- ^^^) = 

Hence, due to cj^aj 7^ 0, the above identity holds if and only if 



a„-„./3/3J - ( — — ) = 0, for j e Ng, |j| = fc, 

/3ez= 



and the claim follows. 



We are now ready to prove the main results of this paper. 

Theorem 2.6. Let k Cz Nq. A non-singular subdivision scheme with symbol 
a{z) and associated parametrization in hi. 12]) with some r G reproduces 
polynomials of degree up to k if and only if 

{D^a){l) = Iml^qjiT) and {D^a){£)=0 for e G S', |j|<fc. 

Proof: The proof is by induction on fc. In the case fc = the claim follows by 
part (z) of Proposition By Proposition 1 1 . 71 it suffices to prove the result for 
the stepwise polynomial reproduction. 

"■<=:" We write any polynomial tt of degree k as Tr{x) ~ c^x^ + Tr{x) 

with TT e llfc-i. Let r > 0. We show that the sequence 

dW = L(tW)= E c,(tM)V^(t«), aez4 
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satisfies S''+^'> = S^S'''> = {7r(tL^+^^), a € Z"}. In fact, due to the induction 
assumption, by (|1.12p and Proposition 12.51 we have 



E -E 



nr+1 



= E c.(t(:+i))'+?(t(:+i)) = .(tri)), ae 

The one but last equality is due to 



a-T _ ,r+l) 



Oi 



,r+l "a 



"=>:" Let j e be such that |j| = k with ji = k for some i ~ 1, . . . , s. Let 
the polynomial ti{x) = and the sequence dW = {7r(tL''^), a e V}. On the 
one hand, by similar arguments as above, we get 

h-j / 1 



h<j 



On the other hand, the definition of tg yields 



h<j ^ 



,r+l 



The polynomial reproduction, i.e. d'-'"^^-' = 5'ad'-''-', implies that 



E 

h<j 



h-j f 1 



a — T 



m 



Thus, by induction for h = j we have 



E/3^ 



'a — mf3 



a T 
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The claim follows from Proposition [^3] ■ 



For convergent schemes we readily get the following result, which is due 
to the fact that for convergent schemes the step-wise polynomial reproduction 
implies polynomial reproduction. 

Corollary 2.7. Let k € Nq. A convergent subdivision scheme with symbol 
a{z) and associated parametrization in hi .12]) with some t ^ W reproduces 
polynomials of degree up to k if 

{D^a){l) ^\m\'q-j{T) and {D^a){e)^{) for £ G S', |j|<fc. 

3 Applications and examples 

It is natural to expect that any shift of the subdivision mask does not effect 
the polynomial reproduction properties of the corresponding scheme, which is 
confirmed by the next result. 

Lemma 3.1. A convergent subdivision scheme Sg. with the symbol a{z) repro- 
duces polynomials up to degree k if and only if so does the shifted scheme 
with the symbol a{z) = z"a{z), a. ^If . 

Proof: By Proposition 12.31 and due to a(l) = m**, we get the following identity 
for the suitable t of and t of 



By Leibnitz differentiation formula and due to the fact that D z°^ , £ G Ng 
evaluated at 1 is equal to qi{(y.) in (|2.1I) we have 



By definition of it is a tensor product polynomial, thus, it suffices to show 
that the following two univariate polynomials in r S M are equal 



T 



T + OL. 




Thus, by Corollary 12.71 to prove the claim we need to show that 




n=0 




a 



n 




The claim follows by the one dimensional result 6, Corollary 5.1]. ■ 
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3.1 Box splines 



An s—variate box spline is given by its symbol 

ae(z)=2^nH^' (3.1) 

eee 

where 9 runs through all the columns of the s x n, rank s matrix O € 2«x" 
with n > s. It is well-known that the subdivision schemes associated with the 
symbols a0{z) are convergent, if the matrix G is such that removing any column 
from Q does not change its rank, see [T] p. 127]. Next results gives the correct 
parametrization for box spline subdivision schemes. 

Lemma 3.2. A subdivision scheme with the symbol 8iQ{z) in (j3.1l) reproduces 
linear polynomials if its associated parameter values are as in (|1.12p with 

Veee eee / 

Proof: The result follows from Corollarv l2.4l and the simple fact that 

D^^aeiz) = r-l-J2 (^J^'^^' n ^ e (C \ {0})^ 



Remark 3.3. In case Q is unimodular, i.e. each s x s submatrix of Q has 
determinant 1 or the integer shifts of the corresponding box splines are 
linear independent and therefore the subdivision scheme associated with slq(z), 
if convergent, is non-singular due to Proposition \1.3\ Hence, the results of 
Proposition \2.3\ and of Theorem \2.6\ hold. 



We consider an example of the 3— directional box splines with 

/ 1 ... 1 ... 1 ... 1 
e= 0... 1... 1 1... 1 

\ k times £ times n times 

and the corresponding symbols ae(z) are denoted by 



In the case k = £ = n = 2, results in [S] imply that the degree of polynomial 
generation is fc = 4. Now, to check the degree of polynomial reproduction we 
first use Lemma 13.21 to compute r that guarantees the reproduction of linear 
polynomials, i.e. t — ^{k + n,t + n) = (2,2). Using Theorem 12.61 with this r 
we see that the scheme does not reproduce polynomials of degree k — 2, since 
gi^i(x) = XI • T2 = 4, but D^^^^') B2,2A^) = 18- 



14 



3.2 Three dimensional example 

Three dimensional examples can also be considered. For example, it is easy to 
show that the subdivision scheme with the mask symbol 



a{zi,Z2,Z3) = 2 



GZ1Z2Z3 



1 + Zl\'^ /1 + Z2\^ f 1 + zA^ /1 + ZiZ2Z3\^ 



5 f 1 + Zl \ /l + Z2\^/l + Z3\^/l + ZlZ2Z3^^ 
1 + Zl\^ f 1 + Z2\ /l + Z3\^/l + 2122^:3 ^ ^ 



— 22 , 

4 V 2 



1 + Zi\^ f 1 + Z2\^ f 1 + Z3\ f 1 + Z1Z2Z3 ^ ^ 



1 + zi\^/1 + Z2\^/1 + Z3Y'^/1 + 2:122^3 



5 

-^ZlZ2Z3 

is convergent. Moreover, using Theorem 12.61 we get that for t — (3,3,3), thus, 
the scheme also reproduces linear polynomials. Since £'(^'°''^^a(l) = 46 7^ 8 • 
92.0,0 (t) = 48, the scheme is not reproducing polynomials of degree 2. 

3.3 Interpolatory schemes 

Exactly the same argument for interpolatory schemes as in [6l Corollary 5.3] 
extends to the multivariate case. Interpolatory schemes are such whose mask 
satisfies 

ao ^ 1, ^mcK ~ O5 ^ , 

and therefore, when convergent, with a limit function interpolating the initial 
data as well as all the data generated through the recursions. 

Let us assume that an interpolatory scheme generates polynomials up to degree 
k. Due to the special structure of the symbol of interpolatory schemes 

a{z) = 1 + ^ ae(z) 

eG£;\{0} 

we get 1^^0(1) = for aU j e No with |j| < k. Then, by Proposition [231 
the suitable choice of t in p.l2p for reproduction of linear polynomials is r = 
0. Corollary 12.71 and the definition of in (|2.1[) imply that the scheme also 
reproduces polynomials up to degree k with this r. Thus, the following result 
holds and confirms that our results reproduce results in [llj . 

Proposition 3.4. A convergent interpolatory scheme reproduces polynomials 
up to degree k if and only if it generates polynomials of degree up to k. 

For the butterfly scheme with the symbol 

a(zi, 22) = 4 • z^^Z2^[7ziZ2B2,2.2 - 2ziBi^3,3 - 222^3,1,3 - 2ziZ2S3,3,i] 
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we have r = (0,0), as expected. Since a{zi,Z2) satisfies sum rules of order 
4 (see, again [5]), the subdivision scheme generates cubic polynomials, it also 
reproduces cubic polynomials by Corollary 12.71 

3.4 v^— subdivision 

The approximating -n/S— subdivision scheme from [16J with the mask symbol 

a{z) — - (^ZlZ2 + Z^^Z^^ + Z^^zl + Z^^Z2 + ZlZ^^ + z^z^^) 
+ ^ (zi^ + Z2+ ZiZ2^) + ^ {Z2^ + Zi+ Z^^Z2) 

satisfies sum rules at most of order 2. The associated dilation matrix M = 

1 2 1 

2 ^ satisfies AP = —3/ and the corresponding refinable function is also 

refinable with respect to the iterated mask a(ziz^^, zfz2^) ■ a{z). By Corollary 
12. 4[ the corresponding scheme reproduces linear polynomials, if r = (0,0). 
Thus, associated refinable function has approximation order 2. 



4 Conclusions 

In this paper we give algebraic conditions on the symbol of a multivariate subdi- 
vision scheme with dilation matrix ml, \m\ > 2, that allow us to determine the 
degree of polynomial reproduction of the scheme. These conditions also yield 
the correct parametrization for any convergent subdivision scheme to guaran- 
tee polynomial reproduction of degree at least 1. This is true in particular for 
subdivision schemes associated with box splines. The restriction of a dilation 
matrix of type ml and tensor product structure of the polynomial in (|2.ip let us 
extend the univariate results in [5] easily to the multivariate setting. We believe 
that this paper is an important first step towards the investigation of polynomial 
reproduction of multivariate subdivision schemes with general dilation matrix, 
which is currently under investigation. 
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